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ABSTRACT 

We present the stellar and gaseous kinematics of an Sb galaxy, NGC 3223, with the aim of determining the vertical and radial 
stellar velocity dispersion as a function of radius, which can help to constrain disk heating theories. Together with the observed 
NIR photometry, the vertical velocity dispersion is also used to determine the stellar mass-to-light (M/L) ratio, typically one of the 
largest uncertainties when deriving the dark matter distribution from the observed rotation curve. We find a vertical-to-radial velocity 
dispersion ratio of cr^lcr^ = 1.21 ± 0.14, significantly higher than expectations from known correlations, and a weakly-constrained 
Ks-band stellar M/L ratio in the range 0.5-1.7, at the high end of (but consistent with) the predictions of stellar population synthesis 
models. Such a weak constraint on the stellar M/L ratio, however, does not allow us to securely determine the dark matter density 
distribution. To achieve this, either a statistical approach or additional data (e.g. integral-field unit) are needed. 
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1. Introduction 


The structure of galaxies, and in particular their mass content 
and distribution, has been extensively investigated in the past 
decades. Studying the mass distribution of galaxies is closely 
related to the disk mass problem: how does one associate the ob¬ 
served stellar light distribution to a certain stellar mass distribu¬ 
tion? Typically this is achieved by scaling the stellar light distri¬ 
bution by a factor (the mass-to-light ratio, M/L) to derive the stel¬ 
lar mass distribution. However, the stellar M/L ratio is not known 
a priori, and a number of routes exist to estimate it (De Jong & 
Bell 2007, Courteau et al. 2014); for instance, via stellar popula¬ 
tion synthesis models (Bell et al. 2001, McGaugh & Schombert 
2014), spiral arms (Fuchs 2003), bar dynamics (Weiner et al. 
2001, Perez et al. 2004), and, most importantly for the present 
paper, using vertical velocity dispersions (Bottema 1997, Kregel 
et al. 2005, Bershady et al. 2010a, Martinsson et al. 2013). These 
methods give only very roughly consistent results, discrepancies 
can run up to factors of a few. A closely related issue is the max¬ 
imum disk hypothesis: can the stellar disk contribution to the 
rotation curve be scaled up to its maximum possible value (van 
Albada et al. 1985)? Using the various methods described above 
to derive the stellar mass, authors have reached contrasting con¬ 
clusions, ranging from close to maximum disk to submaximal 
disks (e.g. Weiner et al. 2001, Martinsson et al. 2013). 

The determination of the stellar mass is important for the 
derivation of the dark matter distribution in galaxies, notably be- 


* Based on observations collected at the European Southern Obser¬ 
vatory, Chile, under proposal 68.B-0588. 


cause the shape of the dark matter density distribution is related 
to that of the stellar mass (de Blok et al. 2008). The dark matter 
distribution has been the topic of a debate in the last years, where 
it has become apparent that there is a discrepancy between the 
dark matter density profiles found in dark matter-only A Cold 
Dark Matter (ACDM) simulations of galaxy formation and the 
observed dark matter density profiles (e.g. de Blok et al. 2001, 
Gentile et al. 2004, van Eymeren et al. 2009, Kuzio de Naray 
& Kaufmann 2011). Recently, various models involving the (un¬ 
known) effect of baryons on the dark matter distribution have 
been proposed, trying to reduce or remove the discrepancy be¬ 
tween simulations and observations (e.g. Oh et al. 2011a, Gover- 
nato et al. 2012, Maccio et al. 2012, Cloet-Osselaer et al. 2012, 
Di Cintio et al. 2014). 

The internal kinematics of galaxies have been used to in¬ 
vestigate a wide variety of topics, including the mass distribu¬ 
tion of galaxies (e.g. Bosma 1981, Persic et al. 1996, Oh et al. 
2011b, Gentile et al. 2013) and galaxy evolution, in particular 
disk heating (e.g. Minchev & Quillen 2006, Gerssen & Shapiro 
Griffin 2012, Sellwood 2014). A possible route to estimate the 
stellar mass is to use the relation between vertical stellar veloc¬ 
ity dispersion and surface density, as done by the DiskMass sur¬ 
vey team (Bershady et al. 2010a). Stellar kinematics can also be 
used to investigate the process called “disk heating”, the increase 
of the random velocity of stars in the galaxy disk with time. By 
measuring the vertical and radial velocity dispersions, cr^ and 
ctr, one can test the validity of theories put forward to explain 
the disk heating phenomenon (Jenkins & Binney 1990, Shapiro 
et al. 2003), e.g. spiral structure (Sellwood & Binney 2002), scat- 
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tering off giant molecular clouds (Hanninen & Flynn 2002), or 
the infall of substructure onto the galaxy (e.g. Kazantzidis et al. 
2009). Note that there is still no consensus concerning the dom¬ 
inant mechanism(s) at play. Gerssen & Shapiro Griffin (2012) 
hnd that the variation of cr^ and ctr with Hubble type suggest 
that spiral structure might be the primary radial heating agent, 
while giant molecular clouds are consistent with being the three- 
dimensional heating agent. These results are based on a handful 
of galaxies, and the values of of cr^ and ctr are derived through 
a multi-parameter fit of the observed velocity dispersion and ro¬ 
tation velocity prohles. In the present paper we use a similar 
approach. 

To tackle simultaneously the issues of disk mass and disk 
heating, as a pilot study we observed the Sb spiral galaxy 
NGC 3223, whose Hubble distance (corrected for Virgocentric 
flow) is about 38 Mpc. NGC 3223 was selected because of its 
southern declination (which made it accessible for observations 
with the VLT), its moderate inclination and relatively smooth 
optical appearance. It was part of the sample of 74 galaxies for 
which Palunas & Williams (2000) made 2-D Ha velocity helds, 
and then derived a rotation curve. We obtained Hi data, to extend 
the rotation curve to larger radii, and stellar kinematics along the 
minor and major axes. NGC 3223 was also part of the sample of 
Kassin et al. (2006a,b) who used the rotation curve of Mathew- 
son et al. (1992) to derive the baryonic and dark matter distribu¬ 
tions in a sample of galaxies. 

The paper is structured as follows: in Section]^ we present 
the data we used and the reduction process, then in Section 
we describe the way we analysed the data. The main results are 
shown in Section |4] and then in Section |5] we draw our conclu¬ 
sions. 

2. Data acquisition and reduction 

2.1. Photometry data 

The BVI frames were taken at La Silla using the NTT telescope 
(New Technology Telescope), mounted with the SUSI-2 detector 
(D’Odorico et al. 1998), during the night of 3'^'^ and 4* of Febru¬ 
ary 2002. Ten frames were taken in each band, with an exposure 
time of 60 seconds (B-band) and 30 seconds (V- and I-band). 
The frames were cleaned of cosmic ray artifacts, de-biased, flat- 
helded and sky subtracted with the MIDAS software package. 
The images were then aligned and combined together. For flux 
calibration, the extinction coefficients and zero points were de¬ 
termined using Landolt standard star helds SA98 and SA104(b). 
The images of these star helds were reduced like the BVI photo¬ 
metric images. Sky subtraction is automatically performed when 
the magnitude is determined in MIDAS. A Hux calibration check 
was done by comparing measured huxes with data found using 
the Aladin sofware package (Bonnarel et al. 2000). Finally, an 
astrometric calibration was done using the Gaia package from 
Starlink. 

The H- and Ks-band images were made with the infrared 
imager-spectrometer SOFI (Son OF Isaac, Moorwood et al. 
1998). Ten frames were taken through the H-hlter and ten frames 
through the Ks-hlter, with an exposure time of 10 seconds. Be¬ 
tween each frame an offset frame was taken to determine the sky 
background. The infrared sky background varies rapidly so, fol¬ 
lowing the SOFI manual, the sky background images were made 
using a jitter template from which a median background is deter¬ 
mined and subtracted from the NGC 3223 frames. Flat-helding 
and illumination correction was done using the illumination and 
coiTection fields for the H- and Ks-band downloaded from the 
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ESO SOFI website page. Flux calibration was performed using 
the standard starhelds SJ9111 and SJ9144 from which the zero 
point and extinction coefficients were determined using a growth 
curve and SNR analysis. The same way as for the BVI frames, 
an astrometric calibration was performed with the Gaia package 
from Starlink. 

2.2. Spectroscopic data 

The spectroscopic observations were executed during the 6*and 
7* of January 2012, for a total exposure time of 6.5 hours. Long 
slit spectra were taken along the major (PA 51 deg) and the mi¬ 
nor axis (-39 deg) of NGC 3223 with the FORS (Focal Reducer 
and low dispersion Spectrograph, Appenzeller et al. 1998) de¬ 
tector mounted on the VLT telescope using a 0.7 arcsec wide 
slit in combination with the Gris_1028z grism. This resulted in 
an instrumental velocity dispersion of ~ 20 km s“'. The covered 
wavelength range is 7663-9335 A. while the seeing fluctuated 
around 0.85 arcsec. Standard data reduction was performed with 
MIDAS. It included bias subtraction, flat-helding, cosmic ray 
removal and wavelength calibration. Finally, the spectra were 
aligned and averaged and sky subtraction was performed. The 
2-D error spectra were computed by taking into account the gain 
and the read out noise of the detector. 

For the data analyses we used ULyS!^(Koleva et al. 2009). 
UlySS is a full spectrum fitting package that fits the observed 
spectrum against a linear combination of single stellar popula¬ 
tion models, parametrised by their ages and metallicities. It takes 
into account the extinction and the imperfections of the flux cal¬ 
ibration (or lack thereof) by including multiplicative polynomial 
to the fit. In this case we used a 6* degree polynomial. As base 
models we use Vazdekis single stellar population models (SSPs) 
(Vazdekis et al. 2010) coupled with the CaT library from Cenarro 
et al. (2001). Those were convolved with the line-of-sight veloc¬ 
ity distribution function (LOSVD). Thus, the free parameters in 
the fit were the age and metallicity of the SSPs, the velocity and 
the velocity dispersion of the LOSVD. We used the wavelength 
range between 8300A and 8680A, which is relatively clear of 
sky lines. 

To access the instrumental broadening and the precision of 
the wavelength calibration we used ULySS to fit Gaussians to 
the reduced lamp calibration image. The instrumental dispersion 
we derived changes from 25 km s“' at 7900A to 15 km s“' at 
9200A. This instrumental dispersion roughly corresponds to the 
22.5 km s ' broadening of the SSP models, so we could directly 
use the output of the program as a measure of the galaxy’s physi¬ 
cal velocity dispersion. The precision of the wavelength calibra¬ 
tion was estimated to be around 1 km s ' . For the radial profile 
analyses we binned the spectra to a minimum mean signal-to- 
noise of 15, except for the outer points where we lowered our 
requirements to a signal-to-noise ratio of 5. 

2.3. Hi data 

The Hi observations were performed at the ATCA (Australia 
Telescope Compact Array) with three different array setups: the 
6km configuration on May 11, 2002, the L5km configuration 
on July 21-22, 2002, and the 750m configuration on August 5- 
6, 2002, for a total of ~ 1700 minutes on source. The correlator 
setup was chosen to have 512 16 kHz channels with a total band¬ 
width of 8 MHz. 


' htttp://ulyss.univ-lyonl.fr 
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Fig. 1. Isophotes of the Ks-band image of NGC 3223 (cyan contours) 
fitted with FitSKIRT (green contours). 


The data were reduced using the Miriad software package 
(Sault et al. 1995). We used standard procedures to flag, cali¬ 
brate, and subtract the continuum. The uv-data were then Fourier 
inverted (using a Gaussian taper with a width in the image plane 
of 10 arcsec and a “robust” parameter of 0.5), and we obtained 
a data cube. The “dirty” cube was Hanning-smoothed, yielding 
a velocity resolution of 6.6 km s“^ and an rms noise of 1.2 mJy 
beam^' (corresponding to a column density of 1 x lO'® atoms 
cm“^). We then CLEANed the data cube and restored it with a 
Gaussian beam of FWHM 36.94 x 23.59 arcsec. 

3. Data analysis 

3.1. Photometry 

The BVIHKs images were further analysed using the software 
packages FitSKIRT (which uses genetic algorithms to fit opti¬ 
cal and NIR images, including radiative transfer calculations, De 
Geyter et al. 2013), and BUDDA (BUlge Disc Decomposition 
Analysis, de Souza et al. 2004, which was designed to perform 
2D bulge-disk decompositions of optical and NIR images). We 
performed a 2D bulge-disk decomposition of the Ks-band im¬ 
age because it traces best the stellar disk mass. We modelled our 
image as a Sersic bulge plus an exponential disk. 

The fitted isophotes with FitSKIRT are shown in Fig. the 
best-fit effective radius of the bulge is 4 ± 0.5 arcsec with a Sersic 
index of n = 1.4 + 0.2, and the best-fit exponential scale-length 
of the disk is 18 + 2 arcsec. We double-checked this value using 
BUDDA and we found a consistent result. The derived apparent 
Ks-band magnitudes of the disk and bulge are 7.85+0.08 and 
9.78+0.08, respectively. The derived inclination angle (averaged 
over the five BVIHKs bands) is 44.1+0.5°. 

3.2. Stellar kinematics 

One of the main goals of the present study is to determine the 
stellar velocity ellipsoid, thus we have to take advantage of our 


observations of stellar kinematics along the major and minor 
axes, mostly following Gerssen & Shapiro Griffin (2012, and 
references therein). Simple geometrical considerations lead to 
the following relation between the observed stellar velocity dis¬ 
persions along the major and minor axes ((Tmaj and cr^in, respec¬ 
tively) and the three components of the velocity ellipsoid: 

o-Lj = o-e^in^'+ (1) 


= 0-Rsin^i + O-zCOS^/ (2) 

where ctr, ctz and erg are the radial, vertical, and tangential 
components of the velocity dispersion, respectively. The inclina¬ 
tion angle is indicated with i. 

These are two equations and there are three unknowns, there¬ 
fore a third equation needs to be added. This is usually done by 
assuming the epicycle approximation, which is valid in the disks 
of spiral galaxies because the orbits of stars are close to circular: 


^ 1 I 

crl 2 i 51nR / 


(3) 


where is the circular velocity of the material in the disk 
at a given galactocentric distance r. Equations [T] and ^repre¬ 
sent a set of three equations and (in principle) three unknowns. 
However, 14 is not necessarily exactly equal to the measured 
tangential velocity Vg because of asymmetric drift: 





sin^/ 



R—Incr 
dR 


2 

R 


1 
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R dVc\ ^dcrl, 
2yc dR) dz 
(4) 


where h is the exponential scale length of the stellar disk 
and is the tilt term, which can be safely neglected in these 
observations of external galaxies (Gerssen et al. 1997, 2000). 

If no gas kinematics are present, the four equations [T] 
and 1^ are used to determine the four unknowns ctr, cr^ and erg, 
and 14 . In our case, we have two gas kinematics observations: 
2D Ha kinematics from Palunas & Williams (2000) and our own 
Hi data (radially more extended but with a lower angular resolu¬ 
tion). We will use the Ha data as an extra piece of information 
and as a consistency check for the circular velocities determined 
from stellar kinematics. 

However, direct inversion of the equations to determine the 
unknowns as a function of radius results in very large uncertain¬ 
ties due to error propagation and quite noisy data to start with. 

Therefore we followed the approach taken by e.g. Gerssen & 
Shapiro Griffin (2012, and references therein), where the radial 
and vertical components of the velocity dispersion are modelled 
as exponentials (in the radial direction): 


o-R = crR,oe (5) 

o-z = (6) 

and the circular velocity as a function of radius (the rotation 
curve) is modelled as a power-law: 


Vc = VoR“ 


(7) 
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Hence the problem reduces to modelling the observed 
o'majCf); gas rotation curve with six param¬ 

eters (ctr o, 0 -^, 0 , Qr, Qz, Vo, and a). 

The central parts of the galaxy are clearly dominated by 
the bulge. Therefore, in our analysis we only consider the radii 
where the disk-to-bulge brightness ratio is larger than 20. Based 
on the disk and bulge parameters derived in section [TT| we thus 
only consider radii larger than 3.2 kpc (17.3 arcsec). 

In a locally isothermal disk the vertical density distribution 
is given by sech^(z/zo) where z is the vertical distance from mid¬ 
plane and zo is the vertical scale-height. In this case, the disk 
surface density 2, and the scale-height zo are linked to the verti¬ 
cal velocity dispersion through the following relation; 


crl - ttGSzo (8) 

As shown by e.g. Bottema (1993), if the disk mass-to-light 
(M/L) ratio is constant with radius, then the equation above can 
be rewritten as follows: 


= nGiMIL)lK,Zo (9) 

and Ik, - is the surface brightness of the stellar disk 

(in the Ks-band) measured here. Assuming the kinematic scale- 
length (for the vertical velocity dispersion) is twice the photo¬ 
metric scale-length, i.e. Uz = 2h (Bottema 1993, Binney & Mer- 
rifield 1998), we can rewrite eq. (|^ to get an independent esti¬ 
mate of the M/L ratio; 


(M/L) = 


cr 


2 

z,0 


nGIoZo 


( 10 ) 


In our analysis, we will investigate models where Uz - 2h, 
Uz - h, and Uz is a free parameter. 


3.3. Hi data cube analysis 

Using the primary beam-corrected data cube, we find a total Hi 
flux of 26.6 Jy km s', corresponding (at a distance of 38.1 Mpc) 
to 9.1 X 10® Mq. The total Hi map is superimposed onto an op¬ 
tical image in Fig.|^ 

Thorough modelling of the data cube of NGC 3223 goes be¬ 
yond the aim of this paper, but we made some simple models to 
have a good understanding of the Hi distribution and kinemat¬ 
ics. We used the code called TiRiFiC (Jozsa et al. 2007) to make 
models of the Hi data cube, which we iteratively compared to the 
observed one (see Gentile et al. 2013 for more details about the 
procedure). The best-fit model is shown in Fig.[^ it gives a very 
good representation of the observations. The data cube. Fig. 
and the total Hi map indicate that NGC 3223 has a regularly ro¬ 
tating Hi disk, see the velocity field shown in Fig.|^ The gaseous 
disk also has an extension to the SE around position (10 21 45, 
-34 17 30), see Fig. This extension is present in the channel 
maps between velocities 2678 and 2704 km s“*. 

Particular attention was given to the derivation of the rota¬ 
tion curve. A first estimate was derived from the tilted-ring fit 
of the intensity-weighted velocity field, then the rotation curve 
was iteratively modified to match the observed data cube, in par¬ 
ticular the position-velocity diagram along the major axis. Each 
data point is assigned an uncertainty that is the maximum be¬ 
tween two quantities: (1) the difference between the two sides 
(receding and approaching), and (2) the velocity resolution cor¬ 
rected for inclination. The best-fit inclination in the inner parts 
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Fig. 2. Total Hi map of NGC 3223 (contours) overlaid onto an opti¬ 
cal DSS image (colour scale). Contour levels are (1.5, 3, 6, 12) xlO^® 
atoms cm“^. Note that the central contour represents a drop in Hi surface 
density. 



Right Ascension (J2000) 


Fig. 3. Hi data cube of NGC 3223. The beam is shown in the bottom 
left of the top left panel (36.9 x 23.6 arcsec). The contours are -3, 3 
(2.5 (t), and 9 mjy beam^'. Negative contours are dashed. The observed 
data cube is shown as thick black contours and the model presented here 
is shown as thin red contours. 
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Fig. 4. Intensity-weighted Hi velocity field of NGC 3223 (contours) 
overlaid onto the total Hi map (colour scale). Contour levels start from 
2704 km s“' in steps of 40 km s“*. 


(52°, but with an uncertainty of about 10° because of the rela¬ 
tively irregular shape of the Hi contours) is slightly higher than 
the one derived from the optical and NIR data (44.1°, see section 
HD- When plotting the two rotation curves together, the Hi data 
points were corrected for the inclination difference; we consider 
the inclination from the optical and NIR data to be more reliable 
because of the less regular shape of the Hi emission. 


4. Results 

The first fit we made was with the six parameters (ctr q, ctz q, ar, 
Uz, Vq, and a) described in Section [T2] fitting simultaneously the 
major axis stellar velocity dispersion profile, the minor axis stel¬ 
lar velocity dispersion profile, the major axis stellar rotation ve¬ 
locity profile, and the circular velocity curve as derived from the 
gas kinematics (Her, the Hi data were used for the mass decom¬ 
position fits but not for the combined velocity dispersion/circular 
velocity fits). The fit quality is very good but some parameters 
are rather weakly constrained. The results of all the fits are tabu¬ 
lated in Table [2 and the fits are plotted in Fig.|^ 

To investigate the real uncertainties on the parameters in a 
better way, we then made a fit where we forced the two kine¬ 
matic scale-lengths (radial and vertical, and respectively) 
to have the same value, as in Shapiro et al. (2003) and Gerssen 
& Shapiro Griffin (2012). The next choice was to impose az-2h 
in order to be able to use eq. ( [T0| l to estimate the stellar M/L ra¬ 
tio. Lastly, we made a fit with only four free parameters (ctro, 
CTz.o, Vo, and a) and fixing a^-h (to have a constant Toomre sta¬ 
bility parameter with radius) and - 2h (like above, to obtain 
an estimate of the stellar M/L ratio). The uncertainties on the fit¬ 
ted parameters were derived by creating 100 randomly generated 
data sets from Gaussian distributions defined by the observed 
data. This means for each realisation we generated a new data 
point, corresponding to each original data point, by randomly 
sampling from the original point with a Gaussian distribution 
with variance equal to our originally computed variance. These 
new data points were assumed to have the same variance as the 
originals. These randomly generated data were then subjected to 
the same fitting procedure (x^ minimisation) as the original data, 
which gives distributions of values for the best-fit parameters and 


hence estimates of their 1-cr uncertainties from Gaussian fits of 
their distribution. 


4.1. Disk heating 

From observations of eight spiral galaxies, Gerssen & Shapiro 
Griffin (2012) find a relation between Hubble type and ctz/ctr 
ratios; it turns out that later type galaxies have lower (Tz/ctr 
ratios. From the relation they derive, the expected value for 
an Sb galaxy like NGC 3223 would be around 0.63 with an 
uncertainty of about 0.2. From the fit with we find 

o'z/fR = 1.21 + 0.14, not consistent with Gerssen’s relation. If 
we force the ratio to be equal to 0.63, we find (Fig.[^ that the fit 
indeed worsens, in particular the major axis velocity dispersion 
and rotation velocity profiles are underestimated. 

It is possible that we overestimated the inclination (44.1°). 
By taking 40.5°, the value of Kassin et al. (2006a), we find a 
best-fit ratio of ctz/ctr = 1.28 + 0.17, still higher than the value 
suggested by Gerssen & Shapiro Griffin (2012). 

It is also possible that the choice of cutoff radius is incorrect. 
Changing the inner cutoff radius to 15 arcsec (instead of 17.3 
arcsec) gives an almost identical value, cTz/ctr = 1.19 + 0.11. 
Changing it to 25 arcsec gives (Tz/ctr = 1.07 + 0.19, slightly 
lower but (as intuitively expected) with a larger uncertainty, be¬ 
cause there are fewer data points to constrain the parameters. 

The result that NGC 3223 is more vertically heated than ex¬ 
pected implies that NGC 3223 is an outlier in the observational 
relation of Gerssen & Shapiro Griffin (2012). The reason for this 
large vertical heating is not obvious, but recent mergers are ex¬ 
cluded because of the regular morphology and kinematics. We 
note that the values we found for NGC 3223 are marginally in¬ 
consistent with the correlation between cr^ and ctz/ctr found by 
Gerssen & Shapiro Griffin (2012). They plot the velocity disper¬ 
sions at half the kinematic scale-length. With o-^ - 98.3 + 7.3 
km s“* and (Tz/ctr = 1.21 + 0.14, NGC 3223 is (marginally) 
inconsistent with the relation found in their Fig. 6. 


4.2. Constraint on the M/L ratio 

The evaluation of constraints on the stellar M/L was done using 
the results of the fits made with the scale length of the vertical 
velocity dispersion distribution fixed at twice the stellar expo¬ 
nential scale length (see Section|^. 

Applying eq. ( [TOj i to the fits using az-2h and az-2aR=2h in 
Table[2gives values of M/L=0.96+0.34 and 1.28+0.46, respec¬ 
tively, which is quite high for the Ks-band. The errors quoted 
here are formal uncertainties based on the fits and error propa¬ 
gation. We used a sech^ vertical distribution with a scale-height 
of 0.9 kpc, based on the Ks-band exponential scale-length and 
using the average between the scale-length-to-scale-height ra¬ 
tios given by Kregel et al. (2002) and Bershady et al. (2010b). 
We used the Ks-band absolute magnitude of the Sun given by 
Pecaut & Mamajek (2013), i.e. 3.30. 

Could other effects contribute to the real uncertainties? One 
possibility is the inclination angle: if we assume 40.5° (see previ¬ 
ous section) the results above become 0.91+0.33 and 1.08+0.39, 
respectively. Changing the inner cutoff radius to 15 arcsec in¬ 
stead of 17.3 arcsec gives (for az-2h and az-2a^-2h) 0.98+0.35 
and 1.25+0.45, thus a negligible variation in the derived stellar 
M/L ratio. Changing the inner cutoff radius to 25 arcsec instead 
of 17.3 arcsec gives 0.75+0.27 and and 1.36+0.49, respectively. 

The expectations from stellar population synthesis range 
from 0.35 to 0.80, using a B-V colour of 0.80 and the relations 
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Table 1. Best-fit parameters of the joint - I 4 fit. The stellar exponential scale-length A is 18 ± 2 arcsec, corresponding to 3.3 ± 

0.4 kpc. 


Fit 

wr.o 
km s“' 

^z,0 

km s“' 

Or 

kpc 

flz 

kpc 

Vo 

km s ' 

a 

All free 

182+29 

140+16 

5.1+L3 

12.7+3.4 

280.6+1.5 

0.064+0.009 

aR=az 

135+12 

162+12 

8.9+0.8 

—Ur 

284.5+1.8 

0.097+0.014 

a^—lh 

129+15 

192+11 

10 .2+1.8 

6.6 (fixed) 

285.9+1.9 

0.102+0.014 

az—2aji—2h 

232+13 

222+8 

3.3 (fixed) 

6.6 (fixed) 

275.8+1.2 

0.053+0.007 


between colour and stellar M/L ratio found in Bell & de Jong 
(2001), Bell et al. (2003), Portinari et al. (2004), Zibetti et al. 
(2009), and Into and Portinari (2013). The B-V colour seems 
to be a good indicator of M/L ratio (McGaugh & Schombert 
2014), even though at NIR wavelengths the dependency of M/L 
on colour is weak. 

Therefore we can conclude that our estimated range of stel¬ 
lar M/L ratios (between 0.5 and 1.7) is consistent with the most 
massive predictions of stellar population synthesis models. We 
have checked that within the outermost point fitted in the ve¬ 
locity dispersion profiles the stellar surface density (even with a 
stellar M/L ratio of 0.5) is much higher (> 10 times) than the Hi 
surface density: thus our results about the M/L ratio are not fal¬ 
sified by the Hi surface density influencing the vertical velocity 
dispersion. 

4.3. Rotation curve fitting 

Using the constraints on the stellar M/L ratio derived in Section 
|4.2| we made rotation curve fits using different dark matter halo 
density distributions. The bulge and disk contributions to the ro¬ 
tation curve were derived from the parameters of the Ks-band 
bulge/disk decomposition of Section [37T| For the bulge we adopt 
the same stellar M/L ratio as for the disk. The Hi disk contribu¬ 
tion was derived from the observed surface density distribution 
(corrected for primordial He), and for the dark matter halo we 
used two different functional forms: 

• the Burkert (1995) halo, known to give good fits to rotation 
curves (e.g. Salucci et al. 2007). The density p{R) as a function 
of radius R has a central core (i.e. a finite central density) and is 
given by: 


PbuAR) 


PoJ^core 

(R + RcoreXR^ + RL.) 


( 11 ) 


where the central density po and the core radius R^ore are the 
two free parameters. 

• the NFW halo (Navarro, Frenk & White 1996), charac¬ 
terised by a central density cusp, whose functional form is the 
following: 


Pnfw(J^) - 


Ps _ 

(R/Rs){l+R/RA^^ 


( 12 ) 


where ps and R^ are the characteristic density and scale-length 
of the dark matter halo. The parameters ps and R^ can also be 
expressed in terms of the virial mass Mvir and the concentration 
c - Ryir/Rs (Rvir is the virial radius). Dark matter simulations of 
structure formation in the Universe show that c and Myir are not 
independent (e.g. Neto et al. 2007). We adopt here the c - Myir 
relation used by Gentile et al. (2007): 


I^vir 


13.6 


My, 


10"Mp 



Article number, page 6 of[^ 


(13) 


The fits are plotted in Fig. the Burkert halo fit is very good, 
and the product poRcore = 232^^3 Mq pc“^ is consistent with the 
value found by Donato et al. (2009) and Gentile et al. (2009) for 
a large sample of galaxies. 

The NFW halo fits are almost as good. In the fit with no con¬ 
straints on the virial mass and concentration the allowed range of 
virial masses spans almost three orders of magnitude. This range 
would have been larger if we had not put a minimum reasonable 
value of 1 for the concentration parameter. The constrained NFW 
halo fit is only slightly worse (jA value ~ 35% higher), we find a 
value of the virial mass of 5.6/^j, x 10'^ Mq, consistent with the 
value of the unconstrained fit. 


5. Conclusions 

We have presented observations of the stars and gas kinematics 
of NGC 3223, an Sb galaxy located at a distance of about 38 
Mpc. 

From the observed stellar velocity dispersion along the ma¬ 
jor and minor axis, from the major-axis stellar rotation veloc¬ 
ity and from the gas kinematics, we determined the vertical and 
radial stellar velocity dispersion as a function of radius. We 
found a vertical-to-radial velocity dispersion ratio of ct^/ctr = 
1.21 + 0.14, significantly higher the relation between (Tz/ctr and 
Hubble type found by e.g. Gerssen & Shapiro Griffin (2012). 

We also used the Ks-band photometry together with the de¬ 
rived vertical velocity dispersion to estimate the mass-to-light 
(M/L) ratio of the stellar disk. We find that the uncertainties are 
quite large, and the allowed range lies between 0.5 and 1.7 (in 
the Ks-band). This lies at the high end of the values predicted 
by stellar population synthesis models. Because of such a weak 
constraint on the stellar M/L ratio we cannot really constrain the 
distribution of dark matter in NGC 3223. However, the method 
proved to yield meaningful results, consistent with other data, 
therefore a statistical approach involving a sample of galaxies 
with similar data and a spread in physical properties will cer¬ 
tainly provide more answers. 
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